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We construct new exact solutions of the focusing Nonlinear Schrodinger equation (NLSE). This 
is a soliton propagating on an unstable condensate. The Kuznetsov and Akhmediev solitons as well 
as the Peregrine instanton are particular cases of this new solution. We discuss applications of this 
new solution to the description of freak (rogue) waves in the ocean and in optical fibers. 
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-Introduction - This research was motivated by inten- 
tion to develop an analytic theory of freak (or rogue) 
waves in ocean and optic fibers. In the recent time 
the simplest and most universal model for description of 
these waves is the focusing NLSE (Nonlinear Schrodinger 
Equation). In application to the theory of ocean waves 
this equation is used since 1968 pj. In nonlinear optics 
it was known even earlier 

The focusing NLSE is the model of the first approxima- 
tion only. For the surface of fluid this model describes the 
essentially weakly nonlinear quasimonochromatic wave 
trains with average steepness not more than 5 • 10 -2 only 
In nonlinear optics its application is also limited to 
the case of waves of small amplitudes (see, for instance 
Q). Nowadays a lot of models specializing the NLSE are 
developed. For the surface waves these are Dysthe equa- 
tions [B|, Q , for the waves in optic fibers are equations 
that include the third time derivatives and more com- 
plex forms of nonlinearity (see for instance 0, HJ). Also, 
freak waves in the ocean were studied by numerical mod- 
eling of exact Euler equations for potential flow with free 
boundary [§, [lfj ■ The behavior of freak waves studied 
by NLSE and by more sophisticated models shows con- 
siderable quantitative difference. Nevertheless, advanced 
improvement of NLSE does not lead to any qualitatively 
new effects. That means that a careful and detailed study 
of NLSE solutions is still very important problem. 

It has been known since 1971 that the focusing NLSE 
is a very special system, integrable by the Inverse Scat- 



tering Method [ll|. Since that time hundreds of arti- 



cles have been devoted to this subject (see for instance 



the monographs [12h14| . In spite of great attention to 
this field, solutions of focusing NLSE in the presence 
of an unstable condensate have not been studied thor- 
oughly enough. However, this case is the most interest- 
ing for developing a theory of freak waves, because there 
is now universal agreement that extreme waves appear 
as a result of modulation instability of quasimonochro- 
matic nonlinear waves. A nonlinear theory describing 
the development of unstable condensate in the frame of 
1-D focusing NLSE has not yet been formulated. More- 
over, even ordinary solitonic solutions propagating over 



unstable condensate have not yet been found. 

Some special cases are already known. In the pres- 
ence of condensate the maximally analytic wave function 
of auxiliary linear Zakharov-Shabat operator has, in the 
right half-plane of the complex variable X — in (where n 
is the spectral parameter), an additional cut on the real 
axis up to the point A = A (here A is the condensate 
amplitude). The soliton adds a simple pole to the wave 
function singularity. Only the cases when the pole is on 



the real axis have been studied. In 1978 Kuznetsov [15 1 
found an important solitonic solution with the pole out- 
side of the cut ReX > A. This solution is a standing 
soliton with a changing amplitude. Later on this solu- 



3 E3]. In 1983 



tion was rediscovered by other authors 
Peregrine [l8| found a solution with the pole exactly at 
the branch point. This is a remarkable instanton ho- 
moclinic solution - the localized events appear from the 
condensate and return back to it. Now there are already 
known multi-instanton solutions (T9L l20l |. The impor- 
tance of these solutions for the development of the theory 
of freak waves is stressed in the article 2l| . Then, the 



solution of Akhmediev with co-authors is widely known 
[22I ]. This solution corresponds to the pole posed inside 
the cut and is periodic in space but local in time. Unlike 
the Peregrine instanton, this solution change the conden- 
sate phase. 

In this article we describe a general solitonic solution 
with the pole at an arbitrary point in the complex up- 
per half-plane of the spectral parameter. This solution 
includes all other known cases and changes the conden- 
sate phase. It is a moving localized soliton with periodic 
nonmonochromatic intrinstic structure. When the pole 
is far away from the cut, but close to real axis, the so- 
lution is similar to Kuznetsov's soliton. However, when 
approaching the cut (but still not on the real axis) the ve- 
locity of soliton propagation and the size of soliton tend 
to infinity, so the Akhmediev limit happens to be very 
singular. 

For constructing this solution we use the advanced 
mathematical technique of the " local d - problem" . This 
method is a modernization of the " dressing method" , de- 
scribed in 00. 
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-NLSE via local matrix d problem - We study solutions 
of the following NLSE 



il Pt - 7;(f xx - (M 2 - A 2 )ip = 



(1) 



with boundary conditions \(p\ 2 — > A 2 at x — > ±oo. Here 
A = A is a real constant. Equation ((T|) is the compat- 
ibility condition for the following ovcrdctcrmined linear 
system for a matrix function "P: 

- = t/*, i_ = (^ + W) * (2) 

Here 



[7 = /A + u 
1 \ /0 p 



1 fM 2 - l^l 2 % 

-M 2 + l^| 2 



2 V 

If ip = 0, system @ has the following solution: 



1 



(3) 



(4) 



Here 



4> = kx + fit, k 2 = X 2 - A 2 
j\ 

Q, = -i\k, q = r : t" 
A + k 

In what follows we assume that k — > A at |A| — > oo. Then 
1 



s v -q-e 
V /T3^2 y-q ■ 



(5) 



Notice that 

fe(-A) = -k(X), q(-X) = -q(X), 0(-A) = 0(A) 
and 

V(-A) = *+(A) (6) 

We consider the <9-problem on the complex A-plane. We 
arc looking for a second order matrix x(A, A, x, t) obeying 
the equation 



Qx 

oX 



(7) 



and normalized by condition x ^ 1 at |A| — >■ oo. Here 
f = %h(X,X)% 1 



/o(A,A) = /+(A,A) 



(8) 



By virtue of ([6]) function / satisfies condition /(A, A) = 
,/ + (A, A). This proves that 

x -(A)= X + (-A) (9) 

The function x has an asymptotic expansion at A — > oo 

(10) 



-, R 
*= 1+ A 



By virtue of © R + = R. The function \ satisfies the 
following system of equations 

dx 



dx 



Ux - XU 



i^ = (XU + W) - X X Uo 



(11) 



Here U and W are given by expressions ([3]). Because 
system (TlTj) is overdeterminatcd, we have the following 
expression for ip: 



tp = A- 2i? (12) 



(12) 



For an arbitrary choice of matrix function /o(A, A) satis- 
fying condition ([5]), function ip is the solution of equation 

-Solitonic solution - Let us choose 

f (\\\- ( «(A,A) 
/o(A ' A) - U(-A,-A) 

The function / now takes the form: 

/(A, A, x, t) = ae^A + a(-X, — A)i?e~ 2 * 

The matrices A, B are degenerate 

A a f3 = 0-0.0(3, B a j3 = c a dp 

a = (l,q), b = (-5,1); c = (q, 1), d=(l,-q) 
and 

("3* • b ) = a a b a = 0, (~cf • "i?) = Cade = 

We now choose 

a(X, A) =CT(A-?7) 

Here <5(A — r/) is the two-dimensional ^-function. Later 
on, the existence of <5-function allows us to work with two 
values of function <f>: 

0(77) = k(r))(x - irjt), <j){-rf) = -~4> 

We now assume that q(if) — q, q{—rj) = — <?• In a neigh- 
borhood of A = rj, X — —rj we expand A and B: 

A = A +A x {X-ri), B = B + B 1 (X + fj) 

* -GO- ^{rX 

A - q f 1 °\ P, - 3 f-1 -2J 
M -—k\-2q Oh Bl --% ( 1 
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Note that now ~~$ = ~$ and d = ar. We will find a 
solution of the 9 - problem in the form of rational 
functions with two poles: 



Note that 



X = l 



V 



(13) 



X — i] A + 77 

where [/, F are constant degenerate matrices: 

U a (i = u a bp, V a p = v a a*fj 

To avoid a singularity we require 

UA = 0, VB Q = 

Substituting (fT3"|) into Q we end up with a simple linear 
system of equations for u a and v a : 

Ua (l + \ce^) - i±M!ce 2 V = a a Ce 2 * 
k rj + rj 



77 + 77 



^-C*e 2 V + (1 + lCe 20 K = 6 Q Ce 



o 20 



By virtue of (15) the solution of equation JT]) is given as 
follows 

ip = A — 2(ui + viq) 

It is convenient to present final answer in terms of a uni- 
formizing variable £: 



Then £ = fle" 



Here 



Now: 



— (xx — -ft) + i — (kx — cut) 



1 A 2 I 

a? = A(R — -) cos(a), 7 = — ^C^ 2 + sin ( 2a ) 

1 A 2 1 

k = A(R+-)sm(a), u = — {R 2 - -^) cos(2a) 

After the proper choice of C we finish with 

Ae 2ia (1 cos(2a) cosh(w + w) - ±(i? 2 + ) cos(u) 



Here 



u = 4> + ' 



cosh(u + w) — i cos(i>) 
. 2 sin(2a) sinh(u + w) + (R 2 - ^y) sin(u) 
cosh(u + w) — i cos(u) 



aerr. — 74, v — 4> — <p* — i(kx — cot) 
1 + R 2 



2i?cos(a) 



, it) = ln(a) 



p -> Ae iia 



\9\ 



A 2 



at x — > —00 
at x — > +00 
at x — >• ±00 



Fig. [T] show a typical solitonic solution at the moment 
of maximum and minimum amplitude. 




FIG. 1: Typical solitonic solution at the moment of maximum 
(blue) and minimum (red) of amplitude. (R — 2, a — j^n) 

In our notation, we obtain Kuztensov's solution in the 
case of real £. When the pole is near the real axis, 
our solution looks like a Kuznetsov soliton moving with 
constant speed. This speed tends to zero in the limit 
— > 0. Thus, there is a broad area of Kuznetsov- 
like solutions near the real £ axis. Fig. [2] shows a typical 
solution corresponding to this case. The blue curve cor- 
responds to the moment of maximum value of |</?| 2 , while 
the red one corresponds to its minimal value. 




(14) 



FIG. 2: Kuznetsov-like solution at the moment of maximum 
(blue) and minimum (red) of amplitude (R — 2,a — j^). 

The Akhmediev case appears for |£| 2 = 1. Near this 
area the solution is a wave train moving with constant 
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speed and without changing its shape. We obtain the 
approximate expression for the envelope s(x,t) of this 
wave train in the case of Arg(^) = j: 



A[ 1 + Q- 



cosh(u + w) 



Here 



(a 2 cosh 2 (u + w) — 4) i 
Q = a(6 + R 4 + 4j - 4a 2 ) 



(15) 



(16) 

For an arbitrary value of Ar<?(£) the envelope s behaves 
similarly but the expression for it is more complicated. 
Fig. EH shows a typical solitonic solution near |£| 2 = 1 and 
its envelope as calculated by (1T5"T) . 




FIG. 3: Soliton and its envelope. (R = 1.01, a = ?) 

The speed and size of the wave train increases without 
bound in the limit case when £ tends to the circle |£| 2 = 1. 
Thus the Akhmediev solution can be understood as a 
special case when the size and speed of the wave train 
tends to infinity. 

Using the presented method one can construct mul- 
tisoliton solutions with different values £i,..,£ n of the 
uniformizing parameter £. By taking all £j — y 1 one can 
obtain " multi-instanton" solutions found recently in ar- 
ticles [illi^. 

-Conclusion - We have discovered new solitonic solu- 
tions of the NLSE. This discovery essentially increases 
the number of candidate solutions for the analytic de- 
scription of rogue waves. We expect that similar solu- 
tions exist in more exact models than NLSE. Since the 
found solutions change the condensate phase they can 
hardly be used for this purpose directly. However, the 
two-soliton solutions with complex conjugated values of 
poles (77,77) are relevant candidates which could compete 
with homoclinic instanton-like solutions. 
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